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Analytic QCD models are those where the QCD running coupling has the physically correct 
analytic behavior, i.e., no Landau singularities in the Euclidean regime. We present a simple analytic 
QCD model in which the discontinuity function of the running coupling at high momentum scales 
is the same as in perturbative QCD (just like in the analytic QCD model of Shirkov and Solovtsov), 
but at low scales it is replaced by a delta function which parametrizes the unknown behavior there. 
We require that the running coupling agree to a high degree with the perturbative coupling at high 
energies, which reduces the number of free parameters of the model from four to one. The remaining 
parameter is fixed by requiring the reproduction of the correct value of the semihadronic tau decay 
ratio. 

PACS numbers: 12.38.Aw,12.38.Cy,13.35.Dx,11.55.Hx 

I. INTRODUCTION 

The running coupling apt(Q^) = ag{Q'^)/TT in perturbative QCD (pQCD) possesses Landau singularities at low 
momenta (0 < < where A^ ~ 0.1 GeV^), where —q^ = and q being the typical momentum of the 
considered process. These singularities are then reflected in perturbative evaluation of spacelike physical observables 
'D{Q^) in terms of powers of a{KQ^) (with k ~ 1). On the other hand, in local field theories with causality [1] the 
observables V^Q^) are analytic functions of in the entire complex Q^-plane, with the exception of the negative 
semiaxis regime < —M^-^^^., where the threshold mass Mthr ~ 10~^ GeV reflects the threshold for production of 
(light) mesons. If I?((5^) is evaluated in terms of powers of a(/i^) (with renormalization scale /i^ = kQ^, k = const. ~ 
1), the aforementioned analytic behavior of the true 'D{Q'^) should be reflected in a{KQ^). In perturbative QCD this 
is not the case, and the Landau singularities of perturbative coupling apt(Q^) should thus be regarded as unphysical. 

On the other hand, analytic QCD models have the running coupling Ai{Q^) (instead of apt{Q^)) which reflects 
more correctly the analytic properties of I?(Q^)- Nonperturbative studies of ghost-gluon vertex and ghost and gluon 
propagators, using Schwinger-Dyson equations [2] and lattice calculations [3|, give QCD coupling a{Q^) with a flnite 
value at = and without Landau singularities at positive Q^. The minimal analytic (MA) model of Shirkov 
and Solovtsov [3HS] is the first and to date the most widely used of the analytic QCD models. It "minimally" 
modifies the perturbative coupling apt(Q^), by eliminating the offending part of the discontinuity function of apt(Q^) 
on the positive axis of Q^-plane (0 < < A'^) and keeping its discontinuity on the nonpositive axis unchanged. 

This results in a new coupling A^^'^\Q'^) analytic in the complex plane with the exception of the nonpositive 
semiaxis — oo < < 0. It turns out that, at high Q^, the MA coupling differs from the perturbative one by terms 
A^^'^'^IQ'^) — apt(Q^) ~ (A^/Q^). Afterwards, other analytic models for Ai{Q'^) have been presented in the literature 
[7HT5]. and they satisfy certain additional constraints at low and/or at high Q^. For further literature and reviews of 
various analytic QCD models, see Refs. [HHIl]- 

Most of the analytic QCD models suffer from one or both of the following problems: 

1. Analytic QCD models with few free parameters usually cannot reproduce the experimental value of the semi- 
hadronic tau decay ratio r,-. For example, MA is a model with this problem, cf. Ref. [3 [17]. The value of r,- is 
at present the most precisely measured low-energy QCD observable, with the squared momentum of the process 

2. The deviation of the analytic coupling Ai{Q'^) from the perturbative coupling apt((5^), at high ^ A^, can 
be appreciable: Ai{Q^) — apt(Q^) (A^/Q^)''' (fc = 1, or 2, ...). This implies that such analytic QCD model 
gives nonperturbative contributions which are at least partly from the ultraviolet (UV) regime ~ {h? /Q"^)^. 
Such UV nonperturbative contributions contravene the Operator product expansion (OPE) philosophy of the 
ITEP (Institute of Theoretical and Experimental Physics) group , which stipulates that the nonperturbative 
contributions come only from the infrared (IR) regime. 
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In Refs. [TH[TH] two analytic QCD models are constructed such that the correct value of Vr is reproduced; however, 
at high the deviations from perturbative QCD are large, namely, A\{Q^) — apt(Q^) ^ (A^/Q^)*^ with fc = 1. 
Further, in Refs. [TH[T3] a method is presented for the construction of higher order analogs of powers a'^^{Q'^) for any 
analytic QCD model. 

In Ref. ^ an analytic coupling AilQ"^) is constructed which comes close to fulfilling the abovementioned ITEP- 
OPE condition, namely, it achieves Ai{Q^) — apt(Q^) ^ {h? /Q"^)^ with fc = 3 at high . However, the higher order 
analogs and the implications for the value of are not investigated in that reference. 

In Refs. [inill^, analytic coupling Ai{Q^) is calculated from certain classes of beta-functions, and such Ai{Q'^) 
fulfills the OPE condition exactly, i.e., \Ai{Q'^) - apt(Q^)| < (A^/Q^)'" for all fc > at high Q^. However, in such 
models it turns out to be difficult to reproduce the correct value of Vr [El [20] . 

In this paper we present a relatively simple analytic QCD model which approximately fulfills the OPE condition 
(i.e., it merges with perturbative QCD to a high degree of accuracy at high \Q'^\ > 10^ GeV^) and, simultaneously, 
reproduces the correct value of r,- . 

In Section [n] we present the model, by motivating it first with a specific simple form of the disconti nuit y function 
Pi{^) = lTaAi{Q^ — —a — ie) in a specific renormalization scheme (/32 — /S^ = ■■ ■ = 0). In Section 



HI 



we fix the 

free parameter by requiring that the model reproduce the correct value of r^. There we also present results of the 
model for the Bjorken polarized sum rule (BjPSR) dBj{Q^) at low and compare them with experimental values. 
In Section [IV| we present conclusions. 

II. MODEL DESCRIPTION 



The perturbative coupling apt(Q ) = ag{Q^)/Ti has singularities (cut) along the real Q axis at Q < Ql {'-^ A^) in 

n 

relation for apt(Q^): 



the complex plane, where Ql is the branching point. Application of the Cauchy theorem gives the following dispersion 



aAQ') = - daf^-^ (1) 

where p^^*^((t) = Im apt(Q^ = —a — ie) is the perturbative discontinuity function, and 77 — +0. The minimal analytic 
(MA) model of Shirkov and Solovtsov [IHS] consists of removing the offending cut at positive (0 < < Q^), i.e., 
for -Ql < cr < 



,(MA),^2^ _ 1 / Pi (^) 



Ar\Q')-zl daJ^^^. (2) 



In analogous way, the higher order couplings An {n = 2, 3, . . .) are constructed in MA analogously ([H El [15] and 
references therein)^ 



A^^^\Q') = - r d.{^^,, (3) 
Jo (f^ + Q ) 

where pi^'^(cr) — Im ap^{Q^ = —a — ie). 

As mentioned before, once the parameter A (the MS A-scale) is adjusted in MA so that the high-energy QCD 
phenomenology is reproduced (A « 400-440 MeV when nf — 3, Ref. [B]), the value of the semihadronic tau decay ratio 
r-r (strangeless and massless) is predicted to be about 0.14, much lower [T31[T7| than the well measured experimental 

value Vr « 0.203 ±0.004. Another possibly unattractive aspect of MA is that its coupling A^^^\Q'^) has singularities 
along the entire nonpositive Q^ axis {Q^ < 0, i.e., ct > 0) in the complex Q^-plane. This does not reflect closely the 
analyticity properties of spacelike observables 'D{Q^) which have nonanalyticity cut along negative Q^ axis starting 
at a negative threshold value —M^^/. Q^ < —M^^^^. (i.e., cr > M^^^^). The value of the threshold mass A/thr is typically 
a (multiple of) mass of light mesons. One possibility to incorporate such behavior in the analytic coupling is to 



The MA couplings A^n^^^ (n > 1) defined here are the MA couplings of Refs. [J] |6l I15| divided by it. 
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eliminate certain IR interval < cr < M^^j^^. of the cut in the dispersive relation ([2]), resulting in a "modified" MA 
(mMA) coupling 



1 



+00 



da 



(4) 



Such type of change was proposed in Refs. [8].^ In Ref. j2T], it was pointed out that the coupling Q is a Stielt- 
jes function, and that, as a consequence, the paradiagonal Fade approximants [M — 1/M]{Q^) of such a coupling 

systematically converge to the exact values A'{^^^\Q^) when the Fade index M increases.^ The latter fact was 
checked numerically, and it was also shown there that the aforementioned Fade approximations are equivalent to 
approximating the mMA discontinuity function by a sum of M delta terms 



- (mMA) 
-Pi 



(^) 



\.f2 \ 



Im apt(Q^ = -cr 



Sn) 



(5) 
(6) 



where s = a / h? ^ Sn — ^^nl ^^-^ /« are positive dimensionless quantities, and Afthr ~ Mi < M2 < • • • . Therefore, 
it was argued that, although a sum of delta functions appears to be a very crude approximation for the (continuous) 
function /3i(cr), it gives increasingly better expressions for the coupling Ai{Q^) when the number M of deltas in the 

sum ^ increases. We only know the approximate behavior of the true pi{a) at high a ^ ~ p'i'^\'^)]i ^^id 

we do not know the behavior in the IR regime a ^ A^. Therefore, as was also argued in Ref. [21] . in the regime of 
low positive a we can expect that parametrization of the true pi (cr) in terms of one or a few delta functions may lead 
to a reasonably realistic behavior of Ai{Q^) at low \Q^\. In this work we use only one delta function in the IR regime 
of ct's: 



pi(cr) = tt/i^A^ S{a - Ml) + e((T - Ml 
= T:fl 5{s - si) + e(s - So) X r 



(pt) 
1 



(^) 



(7) 
(8) 



where s = cr/A^, si = M^/A^, sq = M^/A^, and r'f^\s) = p^^^\a) = Im aptiQ^ ^ -a - ie). We may expect^ 
< Ml < Mq, i.e., the actual threshold mass is Mi. The discontinuity function pi(cr), Eq. ([t]), is depicted in Figs.[l] 
(a), (b), for the choice of two sets of values of the parameters (Mi, Mq and A), which will be motivated in the next 
Section. 
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^ The authors of [Sj used for before eliminating an IR-cut, a "second" MA model 7 , which is obtained by minimal analytization of the 
right-hand side of the following form of the perturbative renormalization group equation (RGE): din a/dln = — jB^a— Pio^ — j320? — ■ ■ ■ ■ 
Their numerical analysis is then performed at the one-loop level. 

^ Pade approximant {N/M]{x) to a function f{x) is, by definition, a ratio Pn(x)/Qm{x) of polynomials of degree N and M, respectively, 
such that the Taylor expansions of [N/M](x) and of f{x) around a; = agree up to (and including) the term ^ a;^+*^. 

^ Although we do not impose the condition < Mi < Mq, it will turn out to be true in our model (see the next Section). 
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FIG. 1: The discontinuity function pi(cr), Eq. ([tJ, for the values of scale parameters A = 0.487 GeV and: (a) Mi = 0.371 GeV and 
Mo = 0.956 GeV (si = Mf/A'^ = 0.581 and sq = M^/K^ = 3.858); (b) Mi = 0.612 GeV and Mo = 1.275 GeV (si = M^/A^ = 1.579 
and So = Mq/K^ = 6.862). For comparison, perturbative discontinuity function Pj'((t) at positive cr is also included as a dotted curve [in 
1S2 = ft = ■ ■ ■ = renormalization scheme (RSch), and with Uf = 3]. 



A physical interpretation of the ansatz ([7|-([8]) is also possible: If the discontinuity function pi{a) is to simulate, 
in the first approximation, the spectral functions pv{cr) — Im — —a — it) of spacelike observables T>{Q'^), the 

delta-term in Eqs. ([7])-([8]) and in Figs.[l]can be regarded as a narrow width approximation (NWA) of the dominant low- 
energy resonance.^ For an application of NWA ansatz directly to the spectral function of the vector current-current 
correlator (<^ Adler function), see Ref. [22^ 

Furthermore, in a recent work of Ref. [25 a similar idea was applied directly to the spectral function wi(cr) of the 
vector current-current correlator. The ansatz for Vi{it) there is similar to the ansatz ([7])-(|8]). However, instead of 
approximating the low-energy regime [a < Mq) by a simple delta function, the known measured low-energy values of 
vi{a) of the ALEPH and OPAL collaborations [23^5] were used. The threshold value of the onset of perturbative 
QCD (Mq), and the perturbative coupling strength A) were then fixed by two conditions: the correct measured 
value of the (strangeless) r^. has to be reproduced, and the higher-twist terms of the resulting Adler function have no 
terms of dimension two (cx l/Q^). The work of Ref. [23] is a generalization and refinement of an earlier work Ref. |27j . 

The dispersion relation applied to the discontinuity function ([7])-(|8]) then gives us the following expression for the 
analytic coupling Ai{Q^): 



+ CX3 



da 



Pi(^) ^ f? 

{a + Q'^) u + si 



Jpt) 



ds 



is) 



(s + u) ' 



where we use the notation 



qVa' 



(9) 



(10) 



In the renormalization scheme (RSch) with the beta function coefficients /32 = /Ss = ■ • • = 0, the discontinuity function 
is explicitly known in terms of Lambert function W±i, Ref. (28j (see also Refs. [521130] on the use of Lambert functions 
in evaluation of the n-loop coupling) 



(s) = Im 



Im 



1 



W-i (s^s-/5o/ci exp[-i7r(/3o/ci - 1)]) 



Cl 



1 + W+i 



{^te''"^"^"' exp[+i7r(/3o/ci - 1)] 



(11) 



where s = cr/A^, and j3o and Ci = j3i/ I3q are the two (universal) coefficients in the renormalization group equation in 
the aforementioned renormalization scheme 



da(Q') 
dlng2 



= -^oa^(l + cia) 



(12) 



specifically /3o = (11 — 2ri//3)/4 and (3i — /3qCi — (102 — 38n//3)/16. It is in this renormalization scheme that we will 
asume the form Q of the discontinuity function, and consequently, the form ^ of Ai{Q^). The A scale appearing 
implicitly in these expressions is the Lambert A, and it is related with the MS scale A at n/ = 3 via 



The coupling apt(Q ) in the /32 = A3 = ' ' 
Ql is the branching point Ql = A^Sb, where Sb = 
in this renormalization scheme is known |28] 

apt(Q') 



A« Aexp(0.32G5) . (13) 
renormalization scheme has a Landau cut along < < Ql^ 'where 



d ^i/*^" a; 0.635 when Uf = 3. The perturbative coupling apt(Q^) 



1 



1 



(14) 



^ For two simpler-minded analytic QCD models involving delta-function in pi{<t), see Refs. 1121 113| . 
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where = \Q'^\ exp{i(j>), the upper signs refer to the case < (j> < +7r, the lower signs to — tt < < 0, and 



z± 



1 

cie 



m 

A2 



exp 



zbvr- 



Cl 



(15) 



It can be numericahy checked that the dispers ion relation (nl holds for apt(Q^) in this renormalization scheme, using 
for the discontinuity function the expressions ( 11 ) at s = cr/A^ > 0, and for —Sh < s = a j < the expression 



r[P'\s) = Im 



1 



Cl 



(16) 



One peculiar feature in most of the analytic QCD models is that at large \Q'^\ (^ A^) the analytic coupling Ai{Q'^) 

3t((3^) by an inverse power of Q^: 

\A^iQ') - apt(Q2)| ^ ( ^) HQ'] » A^) . (17) 



differs from the perturbative coupling apt{Q^) by an inverse power of Q^: 

k 



A2 

In MA [HIS], and in the "second MA" of Refs. [7|, this power is A: = 1, i.e., relatively large deviation. 

However, power deviations (171 are not in accordance with the philosophy of the Operator product expansion 
(OPE) as promoted by the ITEP group |18j . According to that philosophy, all nonperturbative contributions such 
as (A^/Q^)'' to (inclusive) observables originate from the infrared (IR) regimes jQ^j < A^. Consequently, the OPE of 
spacelike inclusive observables is interpreted in this approach as superposition of perturbative contributions coming 
from the ultraviolet regime jQ^j ^ A^ (Wilson coefficients) and nonperturbative contributions coming from the IR 
regime (vacuum expectation values of operators). Once we have power deviations of the type (17), this implies that a 
theory with such Ai{Q^) will give us, in evaluations of inclusive spacelike observables, nonperturbative terms (of the 
type A^/Q^) coming at least partly from the UV regime, thus contravening the ITEP philosophy. 

More specifically, the authors of Ref. argued in the following way that the terms ~ [h?' /Q'^)^ in the deviation 



Eq. (17) contravene the ITEP philosophy of OPE. Namely, let us consider the leading-/3o summation of an inclusive 
spacelike observable 'D{Q^) 



p(LB)(g2) 



dt 
T 



Fv{t)a{tQ^e^) 



(18) 



Here, Fxi{t) is the characteristic function of the observable and C = —5/3. The quantity tQ^e^ is the square of internal 
loop momenta appearing in the resummation. In the UV regime t > 1, the deviation (17) then produces power term 
contributions of UV origin in the observable (see also Ref. [32 ) 



dt „ , ^ 



2\k 



For more discussion on these aspects, we refer to Ref. [20 . The aforementioned OPE condition then implies 



|-4i(Q2)_ap,(Q2)| < 



A! 



(|Q2|»A2;fc = l,2,, 



(19) 



(20) 



This condition can be interpreted as a requirement that the analytic QCD model be formally perturbative (analytic) 
QCD at high momenta. We want to construct here a simple analytic QCD model with the discontinuity function pi{a) 
of the type ([7| which is as close as possible to perturbative QCD (in renormalization scheme ^02 = /33 = • ■ ■ = 0) . 
In such a model, all the measured values of the high-energy QCD observables (with jQ^j > 10"'^ GeV^) are then 
reproduced just as in perturbative QCD. In such a model, we have to ensure only that the values of the well measured 



low-energy observables, particularly r^- (with \Q 



3 GeV ), are reproduced. 



It turns out to be difficult to construct analytic QCD theories which fully respect the condition ( 20 ) , see Refs. [T^l^ . 
We will construct here an analytic QCD model of Eqs. ([7])-(|9|, which only approximately fulfills the condition (20 1. 



We will end up with a simple analytic QCD model where the number of free parameters is reduced to one, and the 
latter parameter will be fixed by imposing the condition of the reproduction of the correct value of r^. We note that 
the model of Eqs. ([7|-([9| at first contains three free dimensionless positive parameters (/f , si, sq), in addition to the 
energy scale parameter A. The parameters /f , si and the scale A will be fixed by requiring: the condition that A; = 3 
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in the relation Eq. (171 instead^ of fc = 1; and the correct perturbative QCD value of ^i(Q^) at = (3mc)^. This 
is approximately the highest value of where the number of active quark flavors can still reasonably safely be kept 
n/ = 3 in perturbative QCD with four- loop RGE- running and three-loop matching conditions [55] J 

A consistent accounting of the threshold effects fully within analytic QCD, when nf i— (nf + 1), is not yet known 
for general analytic QCD models. However, in MA of Refs. ^ i5j, such an accounting can be made systematically 
in the following way (see Refs. [6', and also Refs. [35]): At the positive threshold value Q^j^j. = (kTOc)^, with k ^ 1 
(usually in perturbative QCD: 1 < k < 3), the perturbative coupling can be taken to be continuous to a high degree 
of accuracy: apt{Q^^j.;nf = 3) = aptiQ^Y^^', n f ~ 4)- '^^^ latter condition then relates the scale parameters A{nf = 3) 
and A(n/ = 4). In MA, the discontinuity function is then constructed as p'"^'^\(t) = Im apt(Q^ =^ —a — ie; Uj) with 
71/ ~ 3 for a < {Kiric)'^, and nj ~ A for (kTOc)^ < (t < (kto;,)^. Analogous threshold effects are implemented in pf^^^ 
for the transition n/ = 4 i— > 5 at the threshold txthr = {i^^b)^ ■ The function p"^^\a) is thus step-like discontinuous. 
The resulting (global) MA coupling Ai{Q'^)j constructed via the dispersion relation ([2]) with pf^^ replaced by the 

aforementioned step-wise discontinuous p'^^^\ remains analytic in the entire Euclidean region e C\(— oo,0]. The 
values of A(n/)'s can then be fixed by requiring that the model reproduce the measured values of high-energy QCD 
observables (e.g., with |Q2| > iqI GeV^). 

This threshold matching procedure could be implemented also in the present model (in the renormalization scheme 
/?2 = /^s = ■ ■ ■ — 0) by making pi(cr) accordingly step-like discontinuous at (ithr's. However, by assuming that 
at Qtijj. = (3mc)'^ the present ri/ = 3 analytic QCD model merges with perturbative QCD does not result in any 



appreciable error, due to the condition that in the relation (17) we have fc = 3 [see also the comments later in the 



text, the paragraph just after Eq. (25 1 and the paragraph at Eq. (39)]. 



Throughout this work we assume that we are in the regime of three active quark flavors {nj — 3), and that the 



three quarks m, d and s are (almost) massless. The "approximate perturbative QCD" condition A; = 3 in Eq. (17) can 
be expressed via the following two conditions: 



1 

1 
TT 



dsrris) = ft , (21) 

-Sb 

(pt)/ 



dssr^is) = Sift , (22) 



where Sb = = Q^/A^, Ql being the branch point of apt(Q^) in the complex Q^-plane. Conditions (21|-(22) 

mean that the coefficient of (Q^/A^) and (Q^/A^)^ in the deviation Ai{Q^) — apt{Q^) is zero, respectively. Equations 



(21|-(22| were obtained by subtracting Eq. ^ from Eq. ([T]) 



f2 1 r+so J^^ha) 

apt(Q2) - A(Q') = ^ + - / ds , (23) 

U + Si TT J_ (S + U) 



and expanding in powers of (1/m) = {K^ /Q'^^ 



In addition to the two conditions ([2l|)-(22), there is a condition that the theory merge with the perturbative 
coupling (in the renormalization scheme = — ■ ■ ■ — &i higher renormalization scales p? . Since we do not 
yet know a consistent exact threshold conditions within analytic QCD models, we assume that our analytic QCD 
model coupling Ai{p'^) has the number of active quarks n/ = 3 up to the renormalization scale (RScl) p^ — (3toc)^ 
(with rric w 1.27 GeV) and that at that scale it merges with the value of the perturbative coupling apt((3TOc)^; /32 = 
— ■ ■ ■ — 0;nf — 3) such as implied by the high energy QCD experiments. Specifically, high energy QCD implies 
a(M|,MS) « O.IIO/tt, Rcf. [34]. We then run this value, by perturbative four-loop RGE in MS renormalization 
scheme, down to renormalization scale p^ = (3mc)^, and incorporate quark thresholds at p^ = (3mq)^ {q = b,c) by 
three-loop matching conditions Ref. [55]. We obtain in this way apt = apt((3mc)^, MS, = 3) = 0.07245. Changing 
renormalization scheme from MS to /32 = • • • = by the subtracted form (Ref. 06] , Appendix A there) of the 



® This condition, with A; = 3, is also fulfilled in the model for ^i(Q^) of Ref. [9]. 

^ If we apply the four-loop RGE-running in MS renormalization scheme with three-loop matching conditions according to Ref. [33] at 
= (kuIc)^ and = (Km;,)^, and choose a fixed initial value a(m^, MS; = 3) = 0.12945 at the initial Qq = [which gives: 
a((3m,c)^,MS;n^ = 3) = 0.07245, i.e., the value we use in this paper], we obtain the values a3(M|,MS) = 0.1190,0.1191,0.1193, when 
choosing for the threshold parameter k the values k = 3, 2, 1, respectively. For the quark masses we use the values = 1.27 GeV and 
rub = 4.20 GeV [33]. 



7 



TABLE I: The dimensionless nonnegative parameters si — Mf/A'^ and /i as functions of the cutoff parameter sq = Mq/A^ 
(> 0). The scale A is practically independent of sg: A « 0.487 GeV. The penultimate column shows the leading-/3o resummed 
(LB) contributions to r-r; in parentheses the leading order (LO) contribution. The last column shows the sum of the LB and the 
beyond-the-leading-/?o contribution (LB+bLB) to r^; in parentheses the LO and beyond-leading-order contribution (LO+bLO) 
to r-r (for details on r-r, see the next section). 



so S^ ll r<^^) (r<^°)) ^(LB+bLB) (^(LO+bLO)) 

1.958 0.0000 0.1721 0.2522 (0.1315) 0.2399 (0.1892) 

2.000 0.0121 0.1732 0.2509 (0.1315) 0.2386 (0.1893) 

3.000 0.3117 0.1970 0.2290 (0.1316) 0.2166 (0.1915) 

3.858 0.5812 0.2156 0.2156 (0.1317) 0.2030 (0.1939) 

4.000 0.6267 0.2186 0.2137 (0.1317) 0.2010 (0.1943) 

5.000 0.9523 0.2387 0.2016 (0.1319) 0.1885 (0.1975) 

6.000 1.2861 0.2576 0.1916 (0.1321) 0.1781 (0.2006) 

6.862 1.5788 0.2732 0.1844 (0.1323) 0.1704 (0.2030) 



integrated perturbative QCD RGE (see Ref. [37., Appendix A there) then results in^ 

ai„ = apt((3m,)2;/32 =0,/33 = 0,...;n/ = 3) = 0.07050 



(24) 



As stated above, we require that Ai{fJ,'^) of our analytic QCD model, at /i^ = (3toc)^, agrees with the value Eq. (24 1, 
i.e., the analytic QCD model merges with perturbative QCD starting at the scale = (3mc)^ upwards 



Aiin^ = (3mc)2) = flin (= 0.07050) 



(25) 



One may worry that the replacement of the presented analytic QCD model by perturbative QCD at > (3mc)^ [and 
with = 3 I— > 4 perturbative threshold at = (3mc)^] may not be a good approximation, i.e., that the analytic 
coupling Ai{Q'^) of the theory, say at fixed nj = 3, behaves at high scales Q"^ ~ M| significantly different than the 
perturbative coupling apt{Q^)- It turns out that this is not the case. Namely, if we formally keep fixed n/ = 3 (in order 
not to worry about threshold effects in analytic QCD), the perturbative RGE-running ( in t he renormalization scheme 
^2 = • • • = 0) from the initial value apt((3TOc)2; n/ = 3) = 0.070502 = 0.221487/7r [Eq. ([24])] at the scale ^ (3to^)2 
to the high final scale = M| gives the value apt(M|; n/ = 3) = 0.033694 = 0.105852/7^ The analytic couphng, for 
both representative values of parameter sq used later in this work (sq — 3.858, 6.862), gives the same value 0.221487/7r 
at Q2 = (3mc)^, and almost the same values at = ][.f2. _4^(m|) = 0.105853/7r, 0.105856/7r, respectively.^ This 
strongly indicates that Ai{Q^) of the analytic QCD model presented in this work, at both mentioned values of sq, 
is practically indistinguishable from the perturbative apt(Q^) at scales > {inic)^. This conclusion even gets 
generalized to any higher order couplings of this analytic QCD and the perturbative QCD [see the comments in the 
paragraph at Eq. (39l]. 

Altogether, the three conditions (21), ([22]), and (25) eliminate three of the four otherwise free parameters 
so,si,/i,A^ of our analytic QCD model. We are thus eft with only one free parameter, e.g., the dimensionless 
parameter sg in Eq- ^ for Ai{Q^). In Table [l] we present the numerical dependence of the parameters si,/i on 
Sq. It turns out that the value of the scale A « 0.487 GeV is practically independent of the value of sq, it varies by 
less than 0.1% for the range of the so-values presented in Table [l] the reason being that .4i((3mc)^) = apt((3mc)^) 
(= 0.07050) behaves at such scales practically as perturbative coupling Opt due to "approximate perturbative QCD" 
conditions (21 1-(22 In the last two columns of Table[l]we include various evaluated contributions to the strangeless 
and massless semihadronic tau decay ratio r-r] these aspects will be discussed in the next section. In Table [l] we see 



* We use for the MS beta function /3(a{Q-^)) at = {Srric)^ and nf = 3 the Pade [2/3]^(a) based on the 4-loop polynomial MS beta 

function; if using the latter (polynomial) form, we obtain slightly different value ai^ [ = apt((3mc)^; /32 =0,1^3 = 0, - ■ ■ ,nf = 3)] 
= 0.07054. It appears that, at such relatively large values of apt, the Pade [2/3]^(a) is a better approximation to the full (yet unknown) 

MS beta function. 

^ The values of Ai{Q^), for various values of the parameter sq, differ from the values of the perturbative coupling only at low < 
10 GeV^; for example, the relative difference apt{Q^) / Ai{Q^) — 1 is a monotonously decreasing function of for < 10 GeV^; in 
the interval 1 GeV^ < < 5 GeV^ it falls from 0.02 to 0.0004 when sq = 3.858, and from 0.05 to 0.001 when sq = 6.862. 
For example, for the two input values so = 3.858, 6.862 used later in the text, we have A = 0.48679, 0.48687 GeV, respectively; the 
perturbative value of A (at nj- = 3) is ApQCD = 0.48676 GeV. 

The value of sq = 3.858 results in the reproduction of the central value of the experimental result (A5 = 0, niq = 0)cxp = 0.203 it 0.004 
when the leading-/9o (LB) and the beyond-the-leading-/9o (bLB) contributions are evaluated and added together. 
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that the cutoff parameter sq = Mq/A"^ cannot fall below sq w 1.96 because in such a case si = Mf/A'^ turns out to 
be negative and the coupling acquires a Landau singularity (at — — si > 0). 

III. EVALUATIONS OF INCLUSIVE LOW-ENERGY OBSERVABLES 

In this Section we present the results of evaluation of two inclusive low-energy observables in our discussed model. 
We recall that the high-energy observables (with jQ^j > 10^ GeV^) are reproduced in the model because at such 
energies the coupling practically agrees with the perturbative coupling 

\Ai{Q^)~ap,{Q^)\<{A'/\Q^\f {\Q'\»A'). (26) 

The goal here is to fix the only free parameter sq of the model by requiring reproduction of the measured low energy 
QCD observables. The most precisely measured inclusive low energy QCD observable is the semihadronic tau decay 
ratio Rr, which is the ratio of r(r~ z/T-hadrons(7)) and T{t~ — >■ ^'^-e~^7g(7)). After removing the (well measured) 
strangeness-changing contribution, the color and Cabibbo-Kobayashi-Maskawa factors and the electroweak effects, as 
well as the chirality- violating higher-twist ("quark mass") contributions^^ we obtain the experimental value 

r^iAS = 0, TOg 0)cxp. = 0.203 ± 0.004 . (27) 



For details of this extraction we refer to Appendix B of Ref. |20] and references therein. The quantity rr of Eq. (27) 
is timelike, but it can be obtained from the spacelike massless Adler function dAd\{Q^) — apt(Q^) + C'(apt) by contour 
integration [3B] 



1 

2^ 



+7r 



(1 + e^0)3(l - e"^) dAdi(Q' = rny^) , (28) 



The perturbative expansion of dAdiiQ ) 



oo 



dAdi{Q%t = apt + ^(dAdi)„ap+' • (29) 

n=l 



has been calculated up to C(apt) [5^MT] . i.e., the coefficients ((iAdi)n are known for n = 1,2,3. Here, apt = 
apt(/i^; C2, C3, . . .) is at a chosen renormalization scale fj,^ and renormalization scheme (c2, C3, . . .), where Cj = Pj/Po 
{j ^ 2) are the renormalization scheme parameters. For our purposes it is more convenient to reorganize the expansion 



( 29 1 in terms of the logarithmic derivatives 



«Pt.n+i(A^^) ^ ^J^^^^'t! = <t + 0{a;+') {n = 1,2,...), (30) 



(3Sn\ 9(ln^2^« Pt ' ^ pt 
resulting in "modified perturbation" (mpt) series 



dAdi((9^)mpt = Opt + y^(rfAdi) (31) 

n=l 

The first three coefficients (dAdi)n {n = 1, 2, 3) are known since they can be expressed via (rfAdi)fc's [k — 1, . . . , n). The 
basic idea of evaluation of such leading-twist expressions in general analytic QCD models is to replace (cf. Refs. [1^113) : 
for MA, see also Ref. [H]) 

Opt Ai , apt,„+i i-^ An+i (n = 1, 2, . . .) . (32) 
where An+i are the corresponding logarithmic derivatives in analytic QCD 



^■^ Further, assuming that the chirality-conserving higher twist effects are negligible, i.e., that the gluon condensate (aGG) is approximately 



9 



Therefore, the "modified perturbation" (mpt) series of perturbative QCD (31 1 is replaced in analytic QCD models 
by "modified analytic" (man) series 



-4i 



^(c?Adl)ri-4„4 
n=l 



The known truncated series of Adler function in analytic QCD models is then 

dAd\{Q'^)\^L = Ai + (dAdl)l-42 + • • • {dAdl)N~lA 

with TV = 4, and with the coouplings Ai, A2, ... at renormalization scale 

/i^ = Q2 g^p^^) (C - 1) . 



N 



This series is then inserted into the contour integral ( 28 ) , resulting in the sum 

r[^°+^^°)[Nl =/(A,C)+ ^ (dAdi)„/(^„+i,C) , 

71=1 

where N — 4, and I{An+i,C) are the corresponding contour integrals of Ai = Ai, A2, ■ 



I{An+l,C) = 



1 

2tt 



(l + e»*)3(l-e»*) An+iie^my*) 



(34) 



(35) 



(36) 



(37) 



(38) 



The superscript "LO+bLO" in Eq. (37) indicates that this is a sum of the leading order (LO) term I{Ai,C) and of 



the higher order terms beyond-the- leading-order (bLO). 

The expansions (31 1, (34|-(35) are expansions in nonpower quantities apt,n+i and An+i- The latter are constructed 
by applying logarithmic derivative operators directly on the couplings apt and Ai , and are thus formally linear in the 
latter couplings. This is very convenient also for the application of linear integral transforms on the observables, since 
the linear integral transforms of Ai and An+i usually become simply related and since these transforms respect the 
truncation of the series. In MA of Refs. [HIS], the construction of the higher couplings also has such properties, and 
therefore, the transitions from momentum-transfer (Q^) picture to the energy (s or a) picture and to the distance 
(r) picture become elegant and transparent, especially since the truncation of such series is fully respected by the 
transformations (cf. Refs. [SI [33] ) . 

In this context, we stress that the spacelike observables 'D{Q^), such as Adler function or Bjorken polarized sum 
rule (BjPSR), at higher momentum-transfer scales > (3toc)^ obtain practically the same value when evaluated 
in the perturbative QCD by the truncated (modified) perturbation series d{Q'^)\^^^ [cf. Eq. (31)] or evaluated in the 



presented analytic QCD model by the truncated (modified) analytic series c?((5^)mani Eq. (|35|. This is so because! 



At such the values of apt{Q^) and Ai{Q^) are practically equal, as a consequence of the fulfilled condition 
(17) with fc = 3 [see also the comments in the paragraph just after Eq. (25)]. 



• Applying logarithmic derivatives to the relation (17) with k 
we conclude that 



3, valid in the presented analytic QCD model. 



\An+i{Q^) - apt^n+l{Q^)\ 



A2 



(|g2|»A2;n=l,2,.. 



(39) 



By analogy this implies that at > (3mc)^ the values of apt,„+i (in the renormalization scheme (^2 = 1^3 = ■ ■ ■ — 
0) and An+i{Q^) in the presented analytic QCD are practically equal. For example, if keeping Uf = 3 fixed (in 
order not to worry about the threshold effects), we obtain numerically: a2((3mc)^) = 0.005593 and A2{{3mc)'^) = 
0.005592,0.005587 (when sq = 3.858,6.862), and S2(M|) = 0.0012033 and ^2(M|) = 0.0012033,0.0012034 
(when So = 3.858,6.862); further, asHSmc)^) = 0.000468 and AaiiSnic)^) = 0.000467,0.000464 (when so = 
3.858,6.862), and ^3(^1) = 0.0000442 = Aa^Mj) to the digits displayed (when sq = 3.858,6.862). 

There is yet another, more sophisticated, way of evaluating inclusive spacelike and time-like observables. It is based 
on the knowledge of the leading- ^0 part of coefficients d„ and dn of the inclusive spacelike observable such as dAdi(O^)- 
These leading-/3o (LB) parts can then be summed in any analytic QCD in the form 

dt 



(c^Adi)^^) 



Adl 



{t)Ai{tQ'e'') , 



(40) 
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where C = —5/3. The Euchdean {£) characteristic function F^^^{t) is known [15]. Expansion of expression (jiOl in 
logarithmic derivatives AniQ"^) then reproduces exactly the leading- /3o part of the "modified analytic" expansion (34) 



(41) 



where the expansion of each coefficient ((iAdi)n and {dAdijn of the "mpt" and "pt" expansions, Eqs. (31) and (29 1, in 
powers of Pq is 



nn ^nn 



(42) 



fc=-i 



fc=-i 



In practice, we know the full coefficients (dAdOn for n = 1, 2, 3. Subtracting from them the leading-/3o parts CnlpS then 
gives us a truncated series for the beyond-the-leading-/3o (bLB) contributions. Hence, the LB-resummed expression 
(LB+bLB) for the Adler function, in any analytic QCD model, is 



(^Adl) 



(LB+bLB)/^2N[7V] _ 
man ~ 



dt 



N-1 



Fidiit) AiitQ^e'^) + ^ (rAdi)„A.+i 



where N = 4 and 



(TAdl)„ = (dAdl)„ - cW/?o" 



(43) 



(44) 



When inserting the expression ( 43 1 into the contour integral ( 28 1 , we get the LB-resummed expression for in any 
analytic QCD model 



N-l 



^(LB+bLB),[N] ^ ^(LB) ^ ^ (TAdl)„/(-4„+i , C) , 



(45) 



where I{An+i,C) are given in Eq. (381, and the LB-part is obtained by using the LB-integral (40) in contour integral 

(46) 



(28 1. In Ref. [33] this expression was expressed in terms of the Minkowskian coupling 2li(s) 

d(7 



2li(s) 



in the form 



^(LB)^ r^F^^{t)Mte'ml) 
Jo ^ 



(47) 



where the Minkowskian characteristic function F^{t) was calculated explicitly. Using the relation d2li(s)/(ilns = 
—pi{s)/'K and performing integration by parts, the LB-contribution to (the massless and strangeless) can be 
rewritten in the following form in terms of the discontinuity function Pi((t), which is more convenient for numerical 
evaluations: 



,(LB) _ 



1 



dt 



where 



Fr{t)= / *^F,^(0 



(48) 



(49) 



This expression was used in Refs. [El |20], and explicit expression for the characteristic function Fr{t) is given in 
Ref. 20] (Appendix D). 



Since we use a different normalization, F^(t) here is equal to (t/4) times F^(t) of Ref. |43l . 
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For derivations of and more details on the above identities, we refer to Refs. [T51 [TOl I20j . 

For the evaluations of the (strangeless and massless) r^. in our model, we use the LB-resummed expression ( 45 ) 
and, as an alternative, the more rudimentary expression pT] ), both with = 4. In Table [ij of the previous Section, 
we included the results of these evaluations in the las t tw o columns. We used the renormalization scale param eter 
C = 0, i.e., l/i^l — on the contour integral [cf. Eqs. (36), (38)]. We note that the leading-/3o contribution Eq. (48) 
[and: Eq. (40)] is independent of renormalization scale. From the results of Table|l]we see that, when the value of the 
model parameter sq = Mq/A^ increases, the value of slowly decreases in the LB-resummed approach (LB+bLB) 
Eq. 



(45) 



and slowly increases in the other (LO+bLO) approach Eq. (37). The approach of Eq. (45) reproduces the 
central experimental value ([27| — 0.203 at the parameter value 
0.203 at So = 6.862. 



So 



3.858; the approach of Eq. (37) reproduces 



The running couplings Ai{Q ) and AniQ ) (n ~ 2,3) for positive Q {0 < Q < {3m c) ) are depicted in Figs, 
and [3] for the aforementioned choices of parameter values (sq = 3.858 and 6.862), respectively. In these Figures, the 
perturbative couplings (in the same — Pa = ■ ■ ■ — renormalization scheme) are also presented, for comparison. 
We note that the analytic coupling Ai{Q^) coincides with the perturbative coupling apt((5^) at Q"^ = (3mc)^. The 
conditions (21)-(22| bring the behavior of the analytic couplings very close to those ofperturbative couplings in the 
high momentum regime |g2| > A2 (note: w 0.24 GeV^). This is clearly seen in Figs. [2]and[3] where the perturbative 
couplings Opt, Opt, 2 and apt, 3 (dotted curves) virtually agree with the corresponding analytic couplings Ai, A2 and 
A3 for Q2 down to 2 GeV^. The Landau singularities of the perturbative couplings appear for < Ql where 

the branching point is Ql — '^^/^^ 0.635A^ « 0.1504 GeV^ when nj = 3), and perturbative couplings diverge 
at the branching point — Ql (Landau pole). The discontinuity function pi{a) — lmAi{Q^ = —a — ie) for the 



0.4 



0.3 



0.2 





(a) .vo = 3.858 

- — m\ : 




-^AAiQ") : 






|GeV-J ^ _ ^ V^^^'l 

FIG. 2: The analytic couplings A\{Q^) (full curve), A2{Q'^) (dotted curve) and AziQ'^) (dot-dashed curve) for positive Q-^ {0 < Q'^ < 
(3mc)^) for the choice of parameter sq = 3.858 when: (a) linear scale is used for Q^; (b) logarithmic scale is used for Q^. For comparison, 
the corresponding perturbative couplings apt(Q^), apt,2(Q'^) and apt,3(Q^) are included as dotted curves (in the renormalization scheme 
f^2 = = ■ ■ ■ = 0, with Uf = 3). The couplings A2{Q'^) and apt,2(Q^) were rescaled by factor 4, and the couplings A3{Q^) and apt,3(Q^) 
by factor 16, for better visibility. 



two cases sq = 3.858 and 6.862 is presented in Figs. [T] (a), (b), respectively; it differs very much from the perturbative 
analog Pi*(o') which is also included there. 

A closer look at the two approaches, Eqs. (45) and ([37|), reveals: 



• The simple approach Eq. (37) for Vr in this model requires quite a large IR cut Mq for the MA-discontinuity 
function Afg = A^s^ sa 1.28 GeV which appears to be dangerously close to the mass of the process m-r = 1.777 
GeV; in such a case the scales ~ nir are close to the energy regime < < Mq where the discontinuity 
function pi (a) is parametrized by only one delta function - cf . Fig. [l|b) . On the other hand, the LB-resummed 
approach Eq. (45) requires that the IR cutoff be Mq = A^/sq « 0.96 GeV, roughly half of the mass rrir - 
cf. Fig.jTJa). 



The convergence properties of the truncated "modified analytic" sum Eq. ( 37 ) for r-r show that the last (fourth) 
term is appreciable (ss 0.017). On the other hand, the last (fourth) term in the sum Eq. (45) is significantly 
smaller (« 0.005) - see Table [lT| in the case of the renormalization scale parameter C 



0. 



When varying the parameter C, Eq. (36), away from C — upwards, e.g., in the interval between and 



ln(2) (i.e., on the contour between and 2to^), the result for in the approach of Eq. (37) decreases 



by several percent, while in the leading- /3o resummed approach of Eq. (45) it remains virtually unchanged 



12 



0.5 
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(b) 

*o = 6.862 


\ 4A2(e!) 







0.1 

Q'- iGeV^I 



FIG. 3: Same as in Figs. [2] (a), (b), but now for the choice of the analytic QCD parameter so = 6.862. 



TABLE II: The four terms in truncated analytic expansions (451 and (37 1 for r^, for the values of the sq parameter sq = 3 



and So — 6.862, respectively. The renormalization scale parameter C is varied from ln(3/4) to In 2. For C — 0, the two methods, 
with their respective values of so, reproduce the central experimental value = 0.203. 



method, so 


C 


: LB (LO) 


NLB (NLO) 


N^LB (N^LO) 


N^'LB (N^LO) 


sum (sum) 




ln(0.75) 


0.2156 


0.0018 


-0.0258 


0.0150 


0.2068 


LB+bLB, 





0.2156 


0.0015 


-0.0190 


0.0048 


0.2030 


So = 3.858 


ln(1.3) 


0.2156 


0.0013 


-0.0148 


0.0001 


0.2022 




In 2 


0.2156 


0.0011 


-0.0103 


-0.0032 


0.2031 




ln(0.75) 


0.1458 


0.0229 


0.0238 


0.0145 


0.2070 


LO+bLO, 





0.1323 


0.0308 


0.0224 


0.0175 


0.2030 


So = 6.862 


ln(1.3) 


0.1221 


0.0354 


0.0219 


0.0175 


0.1970 




In 2 


0.1085 


0.0399 


0.0224 


0.0164 


0.1872 



see Table |ll| When moving C to negative values < m^), the two approaches have mutually comparable 

stronger renormalization scale dependence, something to be expected since the model is apparently a simple 
approximation to the true situation for pi{a) at low a < m^. 

The dependence of these results on the renormalization scale /i^ is graphically presented in Fig. |4] We see that 



the leading- /3o-resummed (i.e., LB-|-bLB) evaluation method of Eq. (45), in the depicted renormalization scale range 
0.75m^ < l/^^l < 2to^, gives results for rT-(AiS' = 0, nig — 0) that are significantly less renormalization scale dependent 



than those of the simpler evaluation method of Eq. ( 37 ) 



We can compare these results with the corresponding results in perturbative QCD. We recall that the perturbative 
coupling apt(/i^) and our analytic coupling Ai{^j,'^) come together starting at renormalization scale = (SttIc)'^ 
upwards, and we have apt((3mc)^) — Ai{{3mc)'^) w 0.07050 in the renormalization scheme /32 = f^a = ■ ■ ■ — [this 
value corresponds to the value a(Af|,MS) w 0.119/7r, see Sec. |ll] and to a(m^,/32 = /S^ = ■ ■ ■ ^ 0) = 0.3065/7r]. The 
perturbative QCD evaluation cannot use LB-resummation because Landau poles appear in the LB-integral. Therefore, 
we can compare only with the results of the LO+bLO method in perturbative QCD (using the logarithmic derivatives 
Qpt,n's in the contour integrals) and with the results of the usual perturbative "power" expansion method (using powers 



of Opt in the contour integral). The results for r^. are given in Tabic III for the two aforementioned methods in our 



analytic QCD model, and for the LO-|-bLO and "power" expansion method in perturbative QCD, at various values 



of the renormalization scale parameters C. We see from Table III that the variation \Arr \ when the renormalization 
scale parameter C varies between ln(0.75) and In 2 is for the four methods 46, 198, 84, 103, respectively. The facility 
of unambiguous LB-resummation, which is possible only in analytic QCD models, leads to reduced renormalization 
scale dependence of the result for . On the other hand, in perturbative QCD approaches the use of the logarithmic 
derivatives (pQCD LO+bLO) has the tendency to reduce the renormalization scale dependence of the result for 
in comparison with the use of the power expansion, something already noted in Ref. |44) . We also note that in 



In Ref. [44) . perturbative QCD evaluation of the contour integral of r-r was performed in MS renormalization scheme, and in that 
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0.8 
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1.2 
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1. 



2.0 



Exp[C] 

FIG. 4: The dependence of the evaluation method of Eq. ( |45| (LB+bLB) and of Eq. ( |37| (LO+bLO), for r-r, of the contour renormaUzation 
scale parameter exp(C) = |^^|/mj. Values of the so parameter (3.858 and 6.862, respectively) were adjusted so that, at exp(C) = 1, the 
central experimental value rT-(A5 = 0, niq = 0) = 0.203 is reproduced. 



TABLE III: Values of evaluated by the truncated analytic expansions Eq. (|45| [LB+bLB] and Eq. p?] ) [LO+bLO] (for the 
values of the so parameter sq — 3.858 and 6.862, respectively), as well as values obtained by perturbative evaluations LO+bLO 
(involving Spt.n's in the contour integral) and the truncated "power" expansion (involving powers Opt's in the contour integral), 
in the renormalization scheme /32 = /Ss = • • ■ = 0. At renormaUzation scale (3mc), analytic QCD and perturbative QCD are 
assumed to merge: apt((3mc)^) = Ai{{3mc)^) ~ 0.07050. The results are truncated sums of four terms. The renormalization 
scale parameter C is varied from ln(3/4) to In 2. 



c 


rv: LB+bLB 


LO+bLO 


pQCD LO+bLO 


pQCD power exp. 


ln(0.75) 


0.2068 


0.2070 


0.1893 


0.1856 





0.2030 


0.2030 


0.1873 


0.1828 


ln(1.3) 


0.2022 


0.1970 


0.1850 


0.1801 


In 2. 


0.2031 


0.1872 


0.1809 


0.1753 



perturbative QCD, in order to reproduce the correct value of rr « 0.203, we need a larger value of Opt that would 
correspond to apt(M^,MS) 0.121/7r ([M]); in our presented cases, we have apt(M^,MS) ~ 0.119/7r, and therefore 
perturbative QCD gives too low a value of Tt « 0.18 — 0.19. 

Another inclusive low energy QCD observable is, for example, Bjorken polarized sum rule (BjPSR) dBj(Q^) = 
a{Q^) + pja ^) at low . This is a spacelike quantity whose experimental values are included in the last two lines 



of Table lIV) for three representative values of Q^: 1.01, 1.71 and 2.92 GeV^. These experimental data are based on 
the JLab CLAS EGlb (2006) measurements of the T\~^{Q^) sum rule for spin-dependent proton and neutron 
structure functions gf'" [JS]- The measured quantity T\~^ and the "canonical" BjPRS quantity dej are related with 
each other in the following way: 



rr"(Q2) = da^Bj (5r(xBj,Q')-5'/(a;Bj,Q')) (50) 



9A 



J=2 



where gA = 1.267 ± 0.004 [33] is the triplet axial charge, (1 — dBj) is the nonsinglet leading-twist Wilson coefficient. 



scheme the difference between the renormalization scale dependence of the two approaches LO+bLO (there named: modified CIPT) 
and the power expansion approach (CIPT) was found to be even greater than in the here presented case of the renormalization scheme 
/32 = & = ■ ■ ■ = 0. 
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TABLE IV: Bjorken polarized sum rule (BjPSR) results ciBj(Q^) evaluated with the truncated "modified analytic expansions" 
Eq. (|43| [LB+bLB, with sq = 3.858] and Eq. ([35| [LO+bLO, with so = 3.858, 6.862], for iV = 4 and N = 3. The renormahzation 
scale parameter is taken C — (p^ = exp(C)). The number of active quark flavors is n/ = 3. In brackets the variation of 
the result is given when the parameter C increases from zero to In 2, and when the renormahzation scheme changes (at C = 0) 
from /32 = /33 = ■ • ■ = to MS, respectively. The experimentally measured values are given in the last two lines (see the text 
for details). 



So 


method 




dBj(Q^) : = 1.01 GeV^ 


Q2 = 1.71 GeV^ 


Q2 = 2.92 GeV^ 


3.858 


LB+bLB, iV 


= 4 


0.1973 [+2.1%, -10.0%] 


0.1755 [+7.3%, -7.8%] 


0.1595 [+7.2%, -4.3%] 


LB+bLB, iV 


= 3 


0.2290 [+6.1%, -5.3%] 


0.2066 [+5.5%, -3.1%] 


0.1827 [+4.3%, -0.8%] 


3.858 


LO+bLO, iV 


= 4 


0.2774 [-1.3%, -1.3%] 


0.2234 [-4.2%, -0.6%] 


0.1779 [-4.6%, +0.0%] 


LO+bLO, iV 


= 3 


0.2597 [-8.7%, +1.5%] 


0.2061 [-8.7%, +2.9%] 


0.1650 [-7.5%, +3.3%] 


6.862 


LO+bLO, iV 


= 4 


0.2103 [+9.1%, -0.7%] 


0.1926 [+3.0%, -1.1%] 


0.1668 [-1.0%, -0.7%] 


LO+bLO, iV 


= 3 


0.2184 [+0.0%, -2.2%] 


0.1898 [-4.0%, -0.5%] 


0.1598 [-5.5%, +1.1%] 




exp. 




0.23 + 0.18 


0.13 + 0.12 


0.09 + 0.07 






0.23 + 0.12 + 0.13 


0.13 + 0.09 + 0.08 


0.09 + 0.05 + 0.05 



while "/Q^"' ^ U ^ 2) are the higher-twist contributions. The measured JLab values of with the 

elastic contribution excluded, are [45/. 0.1236 ± 0.0254 for = 1.01 GeV^; 0.1605 ± 0.0195 for = 1.71 GeV^; 
0.1789 ± 0.0112 for = 2.92 GeV^. The values given in the last two Hues of Table |iv] are obtained from these 
values by subtracting from the aforementioned measured values the first higher- twist term /Z4 "/Q^ with the value 
~ Air"(Q = IGeV) = -0.040 ± 0.028 obtained by a three-parameter perturbative QCD fit in Ref. gS]. The 
central value (—0.040) was reconfirmed in Refs. j47,i48, by a fit using the MA [11 [5] approach. In addition, in Ref. [47] 
it was shown, with the perturbative QCD and MA approach, that the the exclusion of the elastic contribution leads 
to strongly suppressed coefficients fi^ '^ at the higher-twist terms ^ with j > 3. In the second line of 



experimental values in Table IV the uncertainties were split into the contribution coming from the uncertainty of the 
measured value of T^'"" and the one from the uncertainty of the fitted value /Lt4~". 

The first two coefficients di and ^2 in the expansion of dej were obtained in Refs. |49| (given there in MS renormal- 
ization scheme); the third coefficient ds is not known exactly, but estimates are known, e.g., Ref. [SD]): in MS scheme 
and at renormahzation scale = Q^, it is ((iBj)3 = 130.; we will use this value. The characteristic function Fjy{t) for 
the leading-/3o resummation for BjPSR was calculated and used in Ref. [12] (on the basis of the known [5T| leading-/3o 
parts of coefficients), and was presented in Ref. [13] ■ This allows us to apply the evaluation methods of Eq. (|43| and 



Eq. (35) in our analytic QCD model (in /32 = /^a = • • • = renormahzation scheme) for dBj(Q^)- The results of t hese 
two methods, at three different low squared momenta (Q^ = 1.01, 1.71, 2.92 GeV^) are presented in Table 



IV 



We 

present in Table IV the results both in the case when the N'^LB (and N'^LO) terms of ~ Aa are included in d-BjiQ^) 



{N = 4 case), and when they are not included {N — 3 case). Further, variations of the results under the change of 
renormahzation scale and scheme are also given. The renormahzation scale was varied from the original fi^ ~ 
(C = 0) to /i^ = 2Q^ (C — In 2). The renormahzation scheme was varied from the original scheme — k = ■ ■ ■ = 
to the MS scheme /Ja = 10.0599 and /Jg = 47.2281 (n/ = 3). 

The change of the renormahzation scheme was performed in the following way (cf. Ref. |13|). The dependence of 
couplings on the renormahzation scheme parameters C2 = 1^2/^0 and C3 = /Sa/Z^o is governed by the partial differential 
equations (pDFs) that are obtained from the corresponding pDFs of perturbative QCD under the analytization rule 



(32) 



a^i(Q^;c2,C3) ^ l^Mi 5 

dc2 2 dx^ U""' dx^ ' ^ ' 

dMQ^;c2,C3) ld^A2 



dc2 2 dx'^ ' 

dAi{Q^;c2,C3) ^ 1 d^Ai 
dc3 12 dx^ 



(53) 
(54) 



where x = /Sq ln(Q^/A^). The left-hand sides of these pDFs are truncated, i.e., terms of ^ (^ A^) are ignored, 
because the truncated series for dB}{Q^) is known only up to ^ A4 (if the aforementioned estimated value of (dBj! 
is used). 



3 



The results in Table IV show that the values of dsjiQ ) obtained with the (LB+bLB) method of Eq. (43) with 
N = 4 are worse than those obtained with = 3, since the renormahzation scale and scheme dependence is in 
general stronger in the A^ = 4 case. This has to do with the numerical behavior of the dBj(Q^) series in the approach 
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TABLE V: The coefficients (dAdi)n and (dBj)n for n — 1,2,3, and their leading-/3o (LB) counterparts, in two renormahzation 
schemes: /Sa = /Ss = • • • = 0, and in MS renormalization scheme. The renormalization scale parameter is C = 0, and n/ =3. 



Quantity 


RSch 


di 


di^«) 


da 


4^^' 


da 


df«) 


Adier 


/32 = /33 = ■ • • = 
MS 


1.640 
1.640 


1.556 
1.556 


7.93 
3.46 


15.71 
15.71 


39.00 
26.38 


24.83 
24.83 


BjPSR 


/32 = • ■ • = 
MS 


3.583 
3.583 


4.5 
4.5 


18.32 
13.84 


32.34 
32.34 


91.13 
52.45 


255.23 
255.23 



TABLE VI: The four terms in truncated analytic expansions (35 1 of dBj(Q^) in our analytic QCD model with so — 6.862,3.858; 
and the terms in the corresponding perturbative expansion in apt,n's ("pQCD LO+bLO"); and the terms in the usual power 
expansion, for = 1.01, 1.71 and 2.92 GeV^. The renormalization scale is chosen as fi^ = Q^; renormalization scheme is: 
102 = A = • ■ ■ = 0; in parentheses, the corresponding values in MS renormalization scheme are given. 



[GeV^l 


method 


dBj(Q^): LO NLO N^LO N^LO 


sum 


1.01 


LO+bLO (so = 6.862) 
LO+bLO (so = 3.858) 
pQCD LO+bLO 
pQCD power exp. 


0.1329 (0.1384) 0.0621 (0.0575) 0.0234 (0.0177) -0.0081 (-0.0047) 
0.1367 (0.1471) 0.0746 (0.0798) 0.0485 (0.0366) 0.0176 (0.0102) 
0.1398 (0.1603) 0.0874 (0.1368) 0.0857 (0.1587) 0.0861 (0.2071) 
0.1398 (0.1603) 0.0700 (0.0920) 0.0674 (0.0832) 0.0658 (0.0858) 


0.2103 (0.2089) 
0.2774 (0.2738) 
0.3989 (0.6628) 
0.3430 (0.4213) 


1.71 


LO+bLO (so = 6.862) 
LO+bLO (so = 3.858) 
pQCD LO+bLO 
pQCD power exp. 


0.1142 (0.1198) 0.0507 (0.0502) 0.0249 (0.0188) 0.0028 (0.0016) 
0.1154 (0.1229) 0.0551 (0.0622) 0.0356 (0.0269) 0.0173 (0.0099) 
0.1162 (0.1266) 0.0584 (0.0769) 0.0454 (0.0596) 0.0360 (0.0506) 
0.1162 (0.1266) 0.0484 (0.0574) 0.0387 (0.0410) 0.0314 (0.0333) 


0.1926 (0.1905) 
0.2234 (0.2219) 
0.2559 (0.3136) 
0.2347 (0.2583) 


2.92 


LO+bLO (so = 6.862) 
LO+bLO (so = 3.858) 
pQCD LO+bLO 
pQCD power exp. 


0.0991 (0.1040) 0.0397 (0.0418) 0.0209 (0.0158) 0.0070 (0.0040) 
0.0994 (0.1045) 0.0411 (0.0474) 0.0245 (0.0186) 0.0128 (0.0074) 
0.0996 (0.1057) 0.0418 (0.0505) 0.0269 (0.0298) 0.0176 (0.0188) 
0.0996 (0.1057) 0.0355 (0.0400) 0.0244 (0.0239) 0.0170 (0.0162) 


0.1668 (0.1656) 
0.1779 (0.1779) 
0.1860 (0.2047) 
0.1765 (0.1857) 



of Eq. (43) in this model, because the fourth term (N'^LB, ^ A4) is comparable or even larger than the third term 
(N^LB, ~ ^3) in this approach. E.g., for ^ 1.71 GeV^ the series is (when C ~ and renormalization scheme 



/3; 



0): dsj (Q^) ~ 0.248 - 0.014 - 0.027 - 0.031 



• . We conclude that the leading-/3o resummed approach of 
Eq. (431 is not working well for dBj{Q^) at low Q^, i.e., that the leading-/?o terms are numerically not representative 
(of at least some) of the perturbative coefficients (<iBj)n {n — 2,3, . . .). The latter fact is shown in Table [v] In that 
Table, the perturbative coefficients (ciAdi)n and their leading-/3o (LB) counterparts are also given; these coefficients 
were relevant in the evaluation of Tt- Note that the leading- /3o coefficients are renormalization scheme independent, 
they depend only on the renormalization scale. Comparing the coefficients (dBj)3 and (dAdOs with their leading-/3o 
counterparts, we can understand why the approach of Eq. (|43| with = 4 is expected to work better in the evaluation 
of Tt than in the evaluation of dBj(Q^). We recall that the N^LB term in this evaluation of a spacelike observable 



V{Q^) is TsAiie'^Q^) where T3 = (d^ 



;7(LB)^ 



cf. Eqs. (43)-(44|. 



On the other hand. Table IV shows that the simpler approach of Eq. (35 1 (LO+bLO), in our model gives results for 
dBj(Q^) that in general get more stable under the renormalization scale and scheme variations when the number of 
terms increases from iV = 3 to = 4, for Q'^ > 1.7 GeV^. The case of very low scale = 1.01 GeV^ is an exception, 
and probably has to do with the fact that our analytic QCD model is not very reliable at such low energies. We note 
that for So = 6.862, where this approach is also used, the threshold masses are relatively high: Mi = -^sIA = 0.612 
GeV and Mq = ^sqA = 1.275 GeV. 

However, in order to see whether we have any better convergence behavior in the evaluation of dBj(Q^) than in 
the perturbative QCD, we should compare with the perturbative QCD evaluation of (iBj(Q^)- The perturbative 
coupling and our analytic coupling merge starting at renormalization scale /.t^ — (3mc)^ upwards, where we have 
apt((3mc)^) = Ai{{3nic)'^) w 0.07050, in the renormalization scheme /32 = /^s = ■ ■ • = (this value corresponding to 
0.119/7r, see Sec.|n| 



In Table 



VI 



we present the values of the evaluated expansion terms for 



the value a(M|, MS) R 

rfBj(Q^) at various Q^, in the approach of Eq. (351 in our model (i.e., in terms of ^„'s), in the analogous approach 
in perturbative QCD (i.e., in terms of apt,„'s), and in the usual power expansion approach in perturbative QCD 
(in powers of apt), using the renormalization scale /i^ = (i.e., C = 0), all in the same renormalization scheme 
P2 — ^3 — ■ ■ ■ — 0- In addition, the corresponding values in the MS renormalization scheme are also given there (in 
parentheses). We see that the evaluated expansions in perturbative QCD for (iBj(Q^) at low momenta ~ 1-2 GeV^ 
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behave much worse than the expansion Eq. (35) in our model; in fact, the third and the fourth terms are roughly the 
same in these perturbative expansions (in apt,n's; and in powers) in the renormalization scheme /32 = /^s = ■ ■ • = 0. 
Consequently, the expansions in perturbative QCD turn out to be very unreliable at such low values of . The 
situation in these perturbative QCD expansions becomes even worse in the MS renormalization scheme; this has to 
do primarily with the fact that the offending Landau cut goes quite far into the positive regime in this scheme; the 
branching point (Landau pole) is at Qf, « 0.388 GeV and 0.627 GeV in the two schemes, respectively. 

The experimentally measured values are (Ref.gS]): 0.23±0.18 for = 1.01 GeV^; 0.13±0.12 for 1.71 GeV^ 



and 



VI 



with 



0.09 ± 0.07 for = 2.92 GeV^; and are given also in Table |IV| Comparing the resuhs of Tables |IV 
these values, we see that the results of both methods, Eqs. (35) and (43), in the presented analytic QCD model, with 
So = 3.858 and sq = 6.862, lie in general above the central experimental values, but in general within the large Icr 
uncertainties of the experimental values. On the other hand, the perturbative QCD results for — 1.01 GeV^ and 

= 1.71 GeV^ are significantly higher than those of the analytic QCD model. Some of the perturbative results lie 
outside the la interval of experimental values, and they show in general significantly worse convergence properties 
than the analytic QCD results. Furthermore, the experimental results indicate the tendency to lower values when 

increases, and this is also the case of all the results in Table IV In the case of evaluation of Bjorken polarized 
sum rule in the present analytic QCD model, the method of Eq. (35), i.e., with no leading- /?o resummation, should 
be regarded as the more reliable one. This is in contrast with the results for r-r where we saw that the method of 
Eq. (45), which involves the leading- /3o resummation, is in our analytic model more reliable and less renormalization 
scale dependent. 



IV. CONCLUSIONS 



We presented a simple analytic QCD model which has initially three dimensionless parameters, and the scale 
parameter A. The model is obtained by parametrizing, in a specific renormalization scheme, the unknown behavior 
of the discontinuity function of the coupling in the IR regime by one delta function. The number of parameters in 
the model reduces from four to one by imposing the requirement that the analytic coupling Ai{Q'^) differ from the 
perturbative coupling apt(Q^) at high only by a small amount ^ (A/Q^)"^. Therefore, the model merges with 
perturbative QCD at high energies 

lO^I > iQi GeV^ to a high de gree of accuracy, and reproduces all the values of 
the high-energy QCD observables (with jQ^I > 10^ GeV^) just like the perturbative QCD does. 

The remaining free dimensionless parameter sq is then adjusted so that the model reproduces the well measured 
strangeless semihadronic tau decay ratio r^- - in the case when the leading- /3o resummation is performed [Eq. (45 1; 



So = 3.858] in the evaluation, and in the case when it is not performed [Eq. ( p37| ); so — 6.862]. The evaluated results for 
the tau decay ratio in the approach of Eq. ( 45 ) turn out to be quite stable under the variation of the renormalization 
scale and show good convergence, not quite so the approach of Eq. (37). On the other hand, the evaluated values 
of the Bjorken polarized sum rule (iBj(Q^) at low momentum transfer < 3 GeV^ in the presented model behave 
reasonably well under the variation of the renormalization scale and scheme, and show a reasonable good convergence, 
if no leading- /3o resummation is performed; i.e., in the case of dBj(Q^) the approach of Eq. (35 1 gives better results 
than the approach of Eq. (43 1. The perturbative QCD evaluations of (iBj(Q^) at such low values of turn out to be 
very unreliable. 

It remains an outstanding problem how to perform in a numerically efficient way the change of the renormalization 
scheme for complex values of squared momentum transfer in analytic QCD models in general, and in the presented 
model in particular. Solution of this problem would shed light on the degree of stability of the evaluated tau decay 
ratio Vr under the renormalization scheme variation. 

Another interesting problem would be to parametrize the unknown behavior of the discontinuity function of the 
coupling in the IR regime by two or more delta functions. This would allow us to fulfill the condition of merging the 
model with perturbative QCD (at |Q^| > 10^ GeV^) to an even higher degree of accuracy ~ (A/Q^)^, and would 
thus allow us to apply and interpret the Operator product expansion technique in such analytic QCD models in an 
analogous way as in perturbative QCD. 
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